Ballistic Transport in Superconducting Weak Links in a Microwave Field 
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Nonequilibrium effects and their impact on a charge trans- 
port in superconducting ballistic weak links biased by an ac 
voltage are investigated within the framework of the Keldysh 
technique. We demonstrate that the microwave field destroys 
the phase coherence during the multiple Andreev reflection 
cycle and leads to the effective cooling of subgap quasiparti- 
cles accelerated due to multiple Andreev reflection. For small 
bias voltages this effect results in a strong supression of both 
the excess current and the conductance of the weak link. In 
the opposite limit of large bias voltages the excess current re- 
mains unaffected. We also demonstrate that a simple Boltz- 
mann kinetic approach becomes inadequate if an ac voltage 
bias is applied to the weak link. 

74.40.+k,74.50.+r,74.80.Fp 

Multiple Andreev reflection leads to excitation of quasipar- 
ticles in voltage biased superconducting weak links Q. As a 
result the quasiparticle distribution function for such systems 
is driven out of equilibrium already at small voltages. In the 
case of a time independent voltage bias charge transport in 
ballistic superconducting weak links in the presence of such 
nonequilibrium effects has been studied by Octavio et al. jl]] 
within the framework of a simple classical Boltzmann kinetic 
equation. This so called OTBK model was then widely used 
in a large number of experimental as well as theoretical works. 

Although the OTBK model provides a transparent physical 
picture of multiple Andreev reflection and dissipative charge 
transport in superconducting weak links it remained unclear 
if (and/or under which conditions) this model is sufficient 
to describe quantum nonequilibrium effects in such systems. 
More recently a rigorous theory of charge transport in ballistic 
superconductor-normal metal-superconductor (SNS) struc- 
tures in the presence of a constant voltage bias has been 
developed by means of the Keldysh technique |||3|. In the 
absence of inelastic relaxation the result for a time-averaged 
dissipative current across the system obtained in [^,^] exactly 
coincides with that derived from the OTBK model thus pro- 
viding a formal justification for the OTBK results Q. 

How far can one go applying the OTBK model to vari- 
ous nonequilibrium effects in superconducting weak links? Is 
the agreement between the results Q and specific for 

ballistic weak links biased by a dc voltage or further general- 
ization of the OTBK model (see e.g. jij]) is possible? 

In this Letter we will study nonequilibrium effects in su- 
perconducting microconstrictions in the presence of a time- 
dependent voltage within the framework of the Keldysh tech- 
nique. We shall determine the distribution function and ob- 
tain the current voltage characteristics (CVC) of the weak 



links. We will demonstrate that photon absorbtion and emis- 
sion processes in the weak link (in combination with multiple 
Andreev reflection) may have a strong impact on the trans- 
port properties of the system leading to effective heating or 
cooling of subgap quasiparticles and to a strong suppression 
of the current in the limit of small bias voltages. The latter 
effect may be of interest for applications of SNS structures 
as microwave detectors. We will furthermore argue that such 
nonequilibrium effects cannot be adequately described within 
the (generalized) OTBK model 0. 

The model. Nonequilibrium effects in inhomogeneous su- 
perconductors are conveniently described by means of quasi- 
classical Greens functions in Keldysh-Nambu space 

M F VG + crsh^-G + h^-Gdi + \K, 6} = 0, (1) 
at at L J 

where G is a 2 x 2 matrix in the Keldysh space consisting of 
retarded, advanced and Keldysh Green functions G R ' A and 
G K . The latter are in turn 2x2 matrices in Nambu space. G 
obeys the normalization condition GoG = — £') where the 
"o" product indicates integration over the internal time vari- 
able. The same integration is also implied in the commutator 
[K,G] where K(r,t,t') = H(r,t)lS(t - t') + it(r,t,t') with 
H(r,t) — i[U(r, t)l — A(r, £)]. Here, U is the scalar potential, 
and A is the off-diagonal pair potential. 

We consider a standard model of a short superconducting 
microconstriction (or bridge): two superconducting bulks are 
connected directly or via a small piece of a normal metal of 
a typical size much smaller than the superconducting coher- 
ence length £o (see e.g. p|). We will assume that the exter- 
nal voltage V(t) applied to the system drops at the junction, 
and the electric field does not penetrate into superconducting 
electrodes. The superconducting phase difference across the 
junction is then defined in a standard way (pit) = 2eV(t)/h. 

Nonequilibrium distribution function. We follow Ref. || 
and construct the Green functions by solving Eqs. (|l|) in each 
superconductor and matching these solutions continuously at 
the contact interface. The quantum kinetic properties of the 
system are completely described by the Keldysh Green func- 
tion G K which can be expressed as G K = G R oh — hoG A with 
h = 1 - 2f - 2/ a 3 |§j. With G K , G R and G A known we have 
constructed / and /. The Fourier transforms of the latter 
with respect to time difference represent the "longitudinal" 
and "transversal" components of the distribution function de- 
scribing respectively energy and charge modes |J, a-j, is the 
Pauli matrix. In equilibrium / is equal to the Fermi distri- 
bution function /o and / = 0. A striking feature of ballistic 
SNS bridges is that inside the N-metal the function / is equal 
to zero even if / is driven out of equilibrium due to multiple 
Andreev reflection in the presence of an externally applied 
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voltage. The key reason for this result is that the electric 
field does not penetrate into the ballistic N-metal (the voltage 
drops only across NS interfaces |^,^]) and therefore does not 
cause charge imbalance responsible for a nonzero / compo- 
nent of the distribution function |g| . As the latter component 
cannot be described by means of a classical Boltzmann equa- 
tion we conclude that the condition / = is an important 
prerequisite for the validity of the OTBK model. 

In the limit of a constant voltage V(t) = V applied to the 
system the distribution function / does not depend on time 
and is given by / = /f (E) = £~ f±U E ) with 

n— 1 

f±U E ) = ]jA(ETleV) [l-A{E^neV)] f {E^neV). (2) 

1=0 

Here A(E) = \ , y R (E)\ 2 is the Andreev reflection probabil- 
ity where (neglecting inelastic scattering) ^y R (E) = (E — 
\/E 2 — A 2 )/ A. "+" and "— " label quasiparticles with mo- 
mentum in and opposite to the direction of the current flow. 
Due to multiple Andreev reflection in the presence of electric 
field "+" quasiparticles are accelerated and the distribution 
function (pj) for energies within the gap increases, i.e. the ef- 
fect of heating of such quasiparticles takes place. It is accom- 
panied by a symmetric effect of cooling of "— " quasiparticles 
which distribution function for subgap energies is suppressed. 

Note that in the absence of inelastic relaxation (which 
causes deviations of A(E) from its BCS value) the expres- 
sion (Q) exactly coincides with that obtained by OTBK jlj. 
This coincidence is by no means surprising: electrons and 
holes suffer no scattering in a clean N-metal and obviously 
can be described by the classical Boltzmann equation. Tak- 
ing into account Andreev reflection at NS interfaces by impos- 
ing proper boundary conditions jl| one arrives at the results 
equivalent to those obtained within the general quantum ki- 
netic analysis ||[|. 

The situation changes in the presence of an additional 
microwave field V(t) = V + Vcosut. In this case elec- 
trons and holes moving in the N-metal can absorb and emit 
photons. Obviously such processes cannot be correctly de- 
scribed by the Boltzmann equation which does not contain 
information about off-diagonal elements of the density ma- 
trix. Here we elaborate the quantum kinetic analysis (^) 
and evaluate the distribution function / which now becomes 
time dependent. In a Fourier representation f±(E,t) — 
Yj7=-oo f±.*( E ) exp(±iKLut) we find 

f±, K (E) = [l-A(E)]fo(E)5 Kl0 (3) 

OO OO 

+ £ m KiK ,(±,k,E) f ±)R ,(E^eV^kHu). 

k — — oo k' — — oo 

Here the terms containing the coefficients 

m K , K > = (±i)" - "' J k Jk-( K - K >)"f R ( E ) h R ( E T K MF (4) 

take care about all possible photon absorption and emission 
processes (the probability amplitude for fe-photon processes 
is given by the Bessel function Jk = Jk(eV /fiuS)). 

Eqs. and (Q) selfconsistently determine the quasiparti- 
cle distribution function for a weak link in the presence of an 
ac voltage. This is one of the main results of the present pa- 
per. It is important to point out that terms with k ^ contain 



information about the phase shift of electron and hole ampli- 
tudes due to photon absorbtion and emission and cannot be 
recovered from the Boltzmann equation analysis Q. 

We have solved Eqs. selfconsistently for various 

values of uj, V, V and T. The results for the time independent 
part of the distribution function /±,«=o are shown in Fig. 1. 
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FIG. 1. Distribution function f±(E) at T = 0.9 T c for a) 
eV — 0, b) eV = 0.1 A and different ui. Solid curves: finite 
microwave field eV/hu — 1, dashed curves: results accord- 
ing to Ref. [4], dashed-dotted curves: zero microwave field, 
dotted curves: Fermi distribution fo- 

Our numerical analysis captures all essential features which 
can be summarized as follows. 

In the case V = (Fig. la)) the microwave field excites 
out of the Fermi surface. At frequencies huj < 2A and small 
microwave amplitudes only single photon processes are im- 
portant and we find an increase of f+(E) — f-{E) = f(E) 
at energies —A < E < and a decrease at < E < +A. 
On top of that f(E) exhibits a peak at |i?|<A. This be- 
havior is a consequence of inelastic quasiparticle scattering 
into (and out of) the bound states in the short constriction, 
E± — ±Acos(ifi/2), due to photon absorbtion. The small 
alternating phase difference ip(t) leads on the one hand to a 
certain smearing of the bound state level around E± — ±A. 
On the other hand photon absorbtion of quasiparticles in the 
E+ state yields an emptying of the latter. Symmetrically, 
quasiparticles with energies E < — A absorb a photon and 
fill up the E- state. It is interesting that due to this reso- 
nance effect the subgap distribution function may be tilt to 
the right which corresponds to an effective cooling. For fre- 
quencies Ttw > 2A the sign of the change in the distribution 
function, A/, is altered and we find an effective heating. In 
this case quasiparticles with E < E- < may be scattered 
into the energy range E > 0. In addition we note here (with- 
out showing a figure) that also for huj < 2A the sign of Af 
may be changed if the microwave amplitude V is larger so 
that multi photon processes become important. 

For nonzero V we observe a nontrivial combination of two 
effects: acceleration of quasiparticles due to multiple Andreev 
reflection (cf. Eq. (Q)) and their excitation by a microwave 
field. The former effect depends on the momentum direction 
whereas the latter is sensitive to the sign of the energy E. As 
a result the structure of the distribution function for subgap 
energies turns out to be quite complicated. This function is 
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shown in Fig. lb) (solid lines) for two different frequencies of 
the microwave field. For small uj (see curve with uj = 0.5A) we 
observe a relative cooling of the "+" branch with respect to 
the case V = (the latter curve is depicted by dashed-dotted 
lines). The physical reason for this effect is clear: photon ab- 
sorbtion and emission processes break the multiple Andreev 
reflection cycle thus preventing from further acceleration of 
quasiparticles. But still the acceleration effect dominates for 
such u> and the "+" branch remains overpopulated. On the 
other hand, for larger values of uj (see curve with uj = 3A) 
photon absorbtion dominates over multiple Andreev reflec- 
tion (for uj > 2A quasiparticles are exited outside the subgap 
energy interval already due to one photon processes) and the 
distribution function becomes closer to that for V = 0. Still 
some dependence of / on the momentum direction remains. 

The above physical picture is correct as long as eV <C A, 
i.e. subgap quasiparticles suffer many Andreev reflections. 
For large voltages eV > 2A only one Andreev reflection is 
possible and the microwave field effect less important. For a 
wide interval of voltages we observe additional subharmonic 
gap structures with the period huj (Fig. lb) or satellite struc- 
tures with SE = ±huj in the vicinity of the usual subharmonic 
gap structures (figures not shown). 

As we already pointed out the problem has been recently 
investigated by Zimmerman and Keck Q within the classical 
Boltzmann kinetic equation approach extending the OTBK 
analysis Q to the ac voltage situation. This approach essen- 
tially deals only with diagonal elements of the density matrix 
and contains no information about off-diagonal phase sensi- 
tive terms. Our results reduce to those of Ref. Q provided 
these terms (m K ^ K r with k / in Eq. (^)) are neglected. 
At low voltages these terms are significant, and the results 
ji| (dashed lines in Fig. 1) considerably (and for huj < 2A 
even qualitatively) deviate from those obtained here. On the 
other hand for larger eV ~ A only a few Andreev reflections 
are possible, off-diagonal elements play a minor role and the 
agreement between our results and those of Ref. Q is better. 

Current-voltage characteristics. The current across the 
weak link is determined by the expression for the Keldysh 
Green function G K in a standard way (see e.g. j^H). For the 
time independent component of this current we find 
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V/R + Iexc(u, V, V) + I S hapiro(uJ, V, V) , 



(5) 



where Rq is the Sharvin resistance of the junction, I exc is 
the additional current due to multiple Andreev reflection and 
I Shapiro = X/fc n ^ fc ,n S(2neV — khoj) represents Shapiro peaks 
at the discrete constant voltages V = (k/2n)huj/e. At low T 
the latter includes subharmonics (n > 1) because the current- 
phase relation strongly deviates from a standard sinusoidal 
form U . Below we shall focus our attention on the microwave 
field effect on the excess current I exc = 2 J dE I K =o(E), where 
similarly to Eq. (^) I K is defined selfconsistently by 

oo oo 

/„(£)= J2 E m — '(+.*. ^) ( 6 ) 

fe — — oo re' = — oo 

x {i {E-eV-khu)5 K , fi + I K ,{E-eV-khuj)} , 

with i (E) = (l/2eR )[A{E) - 1] tanh(£/2fci3T). For V -> 
the Shapiro peaks disappear and Eq. (0) is solved by the 
expression I K {E) = 8 K , In(E) witn @~§1 : 
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Y[A(E-leV) [A(E-neV)-l] tanh[ 



E-neV 
2k B T 1 



(7) 
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FIG. 2. CVC in the limit of constant voltage bias and in 
the presence of an additional microwave field for two different 
frequencies. Shapiro peaks are not shown. 

The I-V curves calculated numerically from Eqs. (JEJ) and 
([]) for different uj are presented in Fig. 2 (Shapiro peaks are 
not shown). For large voltages eV 3> A, hui one observes no 
influence of the microwave field on the excess current. With 
the aid of Eq. (^) one can also demonstrate analytically that 
this result remains valid for any microwave amplitude V. 

The physical reason for this result is transparent. Absorb- 
tion (emission) of a photon in the weak link leads to a shift 
of the energy of quasiparticles by the value hu. If only one 
Andreev reflection takes place (i.e. for eV S> A) this shift is 
unique for all quasiparticles and drops out from the expression 
for the current. The same argument holds for multiphoton 
processes if the total energy shift is smaller than eV. 

In contrast to the above limit, for smaller voltages eV < 2 A 
we find a considerable suppression of the excess current I exc . 
The effect becomes particularly pronounced at large frequen- 
cies uj > 2 A (see Fig. 2). This result is a direct conse- 
quence of microwave induced cooling and heating of sub- 
gap quasiparticles with opposite momenta: the difference 
between "+" and "— " distribution functions decreases (see 
Fig. 1) and so does I exc (since the current I is proportional 
to JdE[f+(E-eV/2)- f-(E + eV/2)}). 

To study this effect further let us calculate the low voltage 
conductance G = 1/Ro +2 J dE G K =o(E) as a function of the 
microwave amplitude. For V — > one can derive from Eq. (n) 
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x {g 0tKl (E-khuj) + G K ,(E-khuj)} , 



(8) 



where g 0jK (E) = -ed[i (E) S K , + I K (V = 0,E)]/dE. 

For V — the solution of Eq. (^) reads 
G K (E) = {l/2Ro)A(E)/[l -A(E)]d[ta,nh(E/2k B Tj[/dES K>0 
and we reproduce the conductance peak G — 1/Ro ~ 
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(W£o)(l/Ro)tanli(A/2fe s T) found in Ref. §. The depen- 
dence G(V) is shown in Fig. 3 for different microwave frequen- 
cies. We find a drastic decrease of G(V) already for small mi- 
crowave amplitudes (eV /2hu) 2 < 1. For larger V the effect 
becomes even more pronounced, multiple Andreev reflection 
turns out to be completely destroyed and the conductance is 
suppressed down to its normal state value 1/Ro- 
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FIG. 3. Zero-bias-conductance G(V) normalized to 
G(V — 0) = 15/i?o versus amplitude V of the applied mi- 
crowave field. Dashed curves: G(V) calculated without phase 
sensitive terms m Kj «/, k/0. 



In this paper we discussed the behavior of voltage biased 
weak links in which case the low voltage "foot" structure of 
CVC is caused by multiple Andreev reflection. Although CVC 
of current biased weak links shows a similar structure, in the 
latter case the "foot" is due to a different physical reason - 
the dc Josephson effect. In many experiments the CVC of 
weak links is measured in the regime intermediate between 
the voltage and current biased limits. Therefore it may be 
quite difficult to judge which of the above physical reasons 
could actually explain experimental results, in particular for 
short superconducting constrictions in which case both the dc 
Josephson current and the low voltage excess current are of 
the same order and have the same temperature dependence. 

The results obtained here suggest that a clear distinction 
between these two mechanisms can be easily made in the pres- 
ence of a microwave field. Indeed in a microwave field (and 
at not very low T) stimulation of the Josephson critical cur- 
rent takes place in superconducting microbridges and SNS 
junctions M whereas the low voltage excess current and the 
conductance become strongly suppressed. We believe that 
these opposite trends can be easily identified experimentally. 

In conclusion, making use of the Keldysh technique we de- 
veloped a microscopic theory of nonequilibrium effects in su- 
perconducting weak links in the presence of an external mi- 
crowave field and demonstrated that these effects may have a 
dramatic impact on charge transport in such structures. 

We acknowledge useful discussions with S.N. Artemenko, 
A. A. Golubov, N.Schopohl and G. Schon. The work was sup- 
ported by SFB 195 of the Deutsche Forschungsgemeischaft. 



Discussion. Our analysis of the microwave-induced effects 
in superconducting weak links has led to a transparent phys- 
ical picture which can be summarized as follows. At low dc 
voltages subgap quasiparticles suffer multiple Andreev reflec- 
tion increasing their energy by the value eV after each traverse 
across the weak link. Quasiparticle states with the momen- 
tum direction in (opposite to) that of the current become 
overpopulated (underpopulated) and the low voltage conduc- 
tance increases Q. The phase coherence is preserved during 
the whole multiple Andreev reflection cycle. In the presence 
of an ac field subgap quasiparticles can also absorb and emit 
photons. These processes destroy the phase coherence within 
a multiple Andreev reflection cycle and prevent quasiparti- 
cles from further acceleration or, equivalently, lead to relative 
cooling (heating) of "+" ("— ") subgap quasiparticles in com- 
parison to the case V = (see Fig. 1). As a result both the 
excess current and the system conductance at low voltages can 
be strongly suppressed (Figs. 2,3). The suppression increases 
with the energy of absorbed and emitted photons fkv and 
the amplitude of the microwave field V in which case many- 
photon processes gain importance. This property makes it 
possible to use SNS structures as spectral microwave detec- 
tors. For larger voltages eV > A subgap quasiparticles are 
accelerated and leave the weak link already after a few An- 
dreev reflection events. In this case the effect of an ac voltage 
becomes less pronounced. In the limit eV S> A only one An- 
dreev reflection takes place and the I-V curve is not sensitive 
to an ac field. 
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